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Abstract: We present a new instrumentation and calibration procedure 
for terahertz time-domain spectroscopic ellipsometry (THz-TDSE) that is 
a newly established characterization technique. The experimental setup is 
capable of providing arbitrary angle of incidence in the range of 15°-85° 
in the reflection geometry, and with no need for realignment. The setup 
is also configurable easily into transmission geometry. For this setup, we 
successfully used hollow core photonic band gap fiber with no pre-chirping 
in order to deliver a femtosecond laser into a THz photoconductive antenna 
detector, which is the first demonstration of this kind. The proposed 
calibration scheme can compensate for the non-ideality of the polarization 
response of the THz photoconductive antenna detector as well as that 
of wire grid polarizers used in the setup. In the calibration scheme, the 
ellipsometric parameters are obtained through a regression algorithm which 
we have adapted from the conventional regression calibration method 
developed for rotating element optical ellipsometers, and used here for the 
first time for THz-TDSE. As a proof-of-principle demonstration, results are 
presented for a high resistivity silicon substrate as well as an opaque Si 
substrate with a high phosphorus concentration. Each sample was character- 
ized by THz-TDSE in reflection geometry with different angle of incidence. 
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1. Introduction 

Terahertz time-domain spectroscopy (THz-TDS) has tremendously grown with a wide range of 
applications fll-El. It is now relatively routine to obtain complex (i.e. real and imaginary) spec- 
tral information, for instance the complex dielectric function, with absolute numerical values 
when performing measurements in a transmission geometry. However, transmission measure- 
ments are not possible on many materials and so far the technique has been difficult to apply 
to many metals, thick or highly doped semiconductors, coatings on thick substrates, substances 
in aqueous solution, and any otherwise opaque compound. Moreover, an outstanding technical 



problem with THz-TDS continues to be the determination of absolute spectral values when 
performing measurements in a reflection geometry |5l-|[6l. A reason why THz time-domain re- 
flection measurements are challenging is that the time-domain technique rests on the ability to 
detect the relative amplitude and phase of a time-dependent electric field of a sample as com- 
pared to a reference. In transmission measurements, transmission through an aperture is used 
as a reference. For reflection based THz time-domain a simple mirror can not easily be used 
as a reference because its surface would need to be positioned within a fraction of a micron in 
exactly the same place as the sample so that the reflected phase is accurate. Such precise in situ 
positioning is challenging. 

Ellipsometry is a well established technique in optical range whereby the measurement of the 
two orthogonal polarization components of light reflected at glancing incidence allows a com- 
plete characterization of a samples's optical properties at a particular frequency |7|. Typically, 
one measures the two orthogonal polarization components by a complete 360° characterization 
of the waves amplitude using rotating polarizers. Importantly, ellipsometry obviates the need 
for measurement against a standard reference sample, and so can provide reliable spectroscopic 
information even when surface morphology is unknown, of marginal quality and/or a reference 
is unavailable. It is also self-referencing, so signal to noise ratios can be very good, as source 
fluctuations are divided out. In order to overcome the technical problems mentioned above 
for THz reflectivity, ellipsometry techniques have been recently revisited for terahertz range 
by using a backward wave oscillator source and a Golay cell power detector [8]. There have 
been a number of attempts to extend ellipsometry to far-infrared frequencies using conventional 
Fourier Transform Spectroscopy technology l9l- fT3ll . Unfortunately the lack of sufficiently in- 
tense sources (in addition to the calibration issues we confront here) has meant that such efforts 
have been challenging, although synchrotron based efforts have made some important contri- 
butions in this regard lfT2l - |[T3ll . Generally these studies are limited to even higher frequencies 
(> 4 THz) than we are here. 

Combining ellipsometry with THz-TDS leads to a new technique called terahertz time- 
domain spectroscopic ellipsometry (THz-TDSE), in which a (sub)picosecond pulse with known 
polarization state is used as a probe to illuminate the sample, and then the modified polariza- 
tion state by the sample is detected upon reflection or transmission. Unlike conventional optical 
ellipsometry, the reflected(transmitted) signal is detected coherently in the time-domain which 
allows one to obtain both amplitude and phase of the light in the two orthogonal directions. 
By transforming the time-domain data into the frequency domain through Fourier analysis, it is 
possible to extract ellipsometric parameter spectra similar to the standard optical spectroscopic 
ellipsometry. However, it should be noted that the instrumentation, signal analysis and calibra- 
tion methods in THz-TDSE would differ from those in the standard optical ellipsometry, and 
all need to be revised accordingly. M. Hangyo and his coworkers, in their pioneering works, 
have demonstrated the potential of THz-TDSE for measuring the complex optical constants of 
a Si wafer with low resistivity, the soft-mode dispersion of SrTiC>3 bulk single crystals and the 
dielectric constants of doped GaAs thin films [14j-[ 16]. There are very few other reports on 
using the THz-TDSE ifTvTl technique, however, there has been recent work in developing high 
resolution THz polarimetry [ 18 1-[ 26 1 for material characterization. None of the proposed THz- 
TDSE experimental setups can provide an easy way of changing the angle of incidence without 
tedious work of optical/terahertz realignments. Moreover, less attention has been paid on the 
compensation of the non-idealities of the optical components and the alignments in THz-TDSE 
through a calibration scheme. 

In this paper, we present a new experimental setup for THz-TDSE in which the incidence 
angle in the reflection mode can change very easily in the range of 15° — 85° with no need 
for any realignment. The same setup can be transformed into transmission mode with the same 



ease. Moreover, we propose a calibration scheme that can compensate for the non-ideality of 
the polarization response of the THz photoconductive antenna detector as well as that of wire 
grid polarizers used in the setup. In our calibration scheme, the ellipsometric parameters are 
obtained through a regression algorithm which we have adapted from the conventional regres- 
sion calibration method developed for rotating element optical ellipsometers, and used here for 
the first time for THz-TDSE. As a demonstration, we present the characterization results for a 
high resistivity silicon substrate and a highly phosphorus doped Si substrate. Each sample was 
characterized with different angle of incidence. 

2. Instrumentation 

Figure [T] illustrates our terahertz time-domain spectroscopic ellipsometry setup. In this system, 
the terahertz components are arranged on two straight arms using optical rails; one arm with 
THz emitter is fixed on the optical table, whereas the other arm with THz detector can rotate 
around a center point where the sample is placed. The sample sits on a second rotation stage. 
Such a configuration provides variable incidence angle (15° < < 85°) in reflection mode, and 
can be easily configurable in transmission mode by aligning the arms along a straight line (Fig. 

Be))- 
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Fig. 1. (a) Schematic of the THz-TDSE setup with variable incidence angle in reflection 
mode that is also configurable in transmission mode (M, mirror; RR, retro-reflector; BS, 
optical beam splitter; OL, optical lens; TL, terahertz lens), (b) Lab setup configured in 
reflection mode at 45° incidence angle and (c) transmission mode. 



The setup uses an 8/ confocal geometry with terahertz lenses made of 
poly-4-methyl-pentacene-l (TPX), which is a terahertz- and optically transparent mate- 
rial. Terahertz lenses have 50.8 mm clear aperture diameter and 100 mm focal length, and 
are less prone to misalignments and polarization distortion as compared to off-axis parabolic 



mirrors. The terahertz beam profile can be approximated as Gaussian. The sample is placed at 
the focal point of the terahertz lenses where an effective flat phase front exist at the beam waist. 
Based on the Gaussian beam assumption and 100 mm focal length for the THz lenses, the f# 
is approximately 9 for 0.2 THz. In the present case, the strong focusing and long wavelengths 
minimize the effects of a spread of incidence angles because an effective flat phase front exists 
at the sample position. This is the opposite approach than has been used in previous infrared 
ellipsometers that were designed to work at high /# [13]. Two identical photoconductive 
dipole antennas with collimating substrate lens are used as THz emitter and detector. Two 
rotatable polarizers are placed in the collimated beams immediately before the detector, and 
after the emitter as shown in Fig. [I] Polarizers were wire grid with wire diameter and spacing 
of 10 jUm and 25 /im, respectively, and field extinction ratio of ~ 40 : 1 at 1 THz. 

The laser source is an 800 nm Ti: sapphire femtosecond laser with pulse duration of < 20 fs 
and 85 MHz repetition rate, which is divided into pump and probe beams. The pump beam 
is guided and focused onto the gap of the emitter photoconductive antenna through free- space 
optics in the usual fashion. However, the probe beam is guided toward the detector antenna on 
the rotating arm by a 1 m long hollow core photonic band-gap fiber (HC-PBGF). Using a fiber- 
coupled THz detector facilitates the rotation of the detector arm with no need for realignment of 
the optics for each incidence angle. The advantage of HC-PBGF is that it exhibits extremely low 
nonlinearity, high breakdown threshold, zero dispersion at the design wavelength, and negligi- 
ble interface reflection (27). Near zero dispersion of such fibers elliminates the need for optical 
pulse pre-chirping, which makes the optical setup and the alignments considerably simpler than 
using convensional fibers with pre-chirping. 

The fiber used in our setup is commercially available HC- 800-2 from NKT Photonics. This 
fiber, fabricated by a stack- and-draw method, has a low-loss transmission band-gap of 100 nm 
centered around 820 nm. The core diameter is 7.5 jUm, and the mode-field diameter is 5.5 fim. 
Attenuation through the low-loss region is reported as < 250 dB/km. The fiber is connected 
from one side to a fiber launcher with a focusing lens and a 3 -axis positioner, and from the 
other side is connected to a fiber-coupled focusing optics that is attached to the THz detector 
compartment. 

The temporal terahertz pulse is acquired by a lock-in amplifier during a time window within 
which the time delay between terahertz pulse and the sampling probe laser is swept by continu- 
ous movement of the retro-reflector. The pump beam is mechanically chopped, and the lock-in 
amplifier is synchronized to the chopping frequency. The acquired temporal signal is then taken 
into the frequency domain through a Fourier transform. Fig. [2] shows a typical THz pulse de- 
tected in transmission mode without any sample. In Fig. [2] one can compare the time-domain 
pulse shape and its spectrum detected when the probe beam is delivered to the detector by a 1 m 
long hollow core photonic band-gap fiber vs. free-space coupling. Fig. [3] shows the THz pulse 
and its corresponding spectrum in reflection using fiber coupling for various incidence angles 
without any optics realignment after changing the incidence angle. One can see that in this 
setup, any changes in the overall intensity at different angles are very small. In the reflection 
measurements a silver mirror was used as a reflector. 

3. Formulation 

Real- valued ellipsometric parameter spectra A) are obtained as a result of the spectroscopic 
ellipsometry measurement. These parameters are related to the complex ratio of the polarization 
eigenstates corresponding to the p and s components of the electric field vector (p- and s-axis 
are defined perpendicular to the propagation direction, and parallel and perpendicular to the 
incidence plane, respectively). Specifically, tan 1 ! 1 and A are defined as the absolute value and 
the phase of the complex ratio, respectively, when the incident polarization state is linear and 
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Fig. 2. Typical THz pulse detected in transmission mode when the probe beam is guided 
through (a) 1-m-long hollow core photonic band-gap fiber, (b) free-space optics. Insets 
show the corresponding pulse spectra. 
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Fig. 3. (a) Time-domain THz pulse and (b) its corresponding spectrum in reflection mode 
for various incidence angles. A silver mirror was used as the reflector. 



at 45° azimuthal angle |7|. Once *P and A spectra are obtained from experiments, the optical 
constants and layer thicknesses can be generally extracted by fitting proper dielectric functions 
and optical models to the measured \F and A parameters. In a specific case of a bulk sample the 
complex refractive index (n — jk) can be obtained directly from the ellipsometric parameters 
using the expressions l28l 



a 2 -k 2 = sin 2 6 



Ink = sin 2 6 



tan 2 (cos 2 2<£ - sin 2 2<£ sin 2 A) 
(1 + sin2<£cosA) 2 
tan 2 0sin4<E> sin A 



(1) 



(l + sin2<£cosA) 2 ' 



where is the incidence angle, and tan<£ =1/ tan^. 

Figure [4] illustrates a schematic of the optical configuration of the THz-TDSE in reflec- 
tion mode. A biased photoconductive dipole antenna (Emitter) generates the incident field 
Ef = [Ei p Ei s ] with close to linear polarization state and azimuthal angle around 45°. After pass- 
ing through the polarizer, the polarization state becomes highly linear at the azimuthal angle of 
the polarizer <j> p . The reflected polarization state is then modified according to the ellipsometric 
parameters (*P, A) of the sample. The p- and ^-component of the reflected polarization state are 
detected by setting the azimuthal angle of the analyzer at (j) a = 0° and 90°, respectively. 
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Fig. 4. Schematic of the optical configuration of THz-TDSE in reflection mode. 

In frequency domain, the Jones vector of the light beam at the THz detector, E r = [E rp E rs ] as 
shown in Fig.[4j is determined by the product of the Jones matrices of the polarizer J^, sample 
Js, and analyzer J a as Q 



E r = R(-0 fl ) J A R(0 fl )J s R(-0 p ) J P R(0 p )E/, 



(2) 



where R(0) is the rotation matrix with angle (f). The Jones matrices for the sample and the 
polarizer/analyzer are given as Q 



sin( x I / )exp(jA) 

cos (¥) 



Jp(A) 



1 

7] P(A) exp(j5 P(A) ) 



(3) 



where t\p(a) is the extinction ratio of the polarizer (analyzer), defined as the ratio of the trans- 
mitted electric field when the polarizer transmission axis is perpendicular and parallel to the 
electric field direction, and 8pu) is the phase retardance. For an ideal polarizer/analyzer the 
extinction ratio is zero (T]pu) = 0). 



In practice, the response of the THz detector is far from ideally linear polarization due to 
the non-ideality of the antenna geometry, and optical and terahertz alignments (29][T9]]. There- 
fore, the polarization response of the THz detector should be taken into account by introducing 
a polarization state vector Pj = [P p P s ], which is generally frequency dependent. In the fre- 
quency domain, the spectrum of the output signal of the photoconductive antenna detector can 
be expressed by an inner product as l30l 

I = F d (co)-E n (4) 

where / is the Fourier transform (spectrum) of the THz detector output signal, E r is the electric 
field vector impinging on the detector as shown in Fig. |4| and CO is the angular frequency. We 
define a complex parameter p as the ratio of the detector output spectrum when the transmission 
axis of the analyzer is along p- and s-axis which corresponds to (j) a = 0° and 90°, respectively, 
as 

At/ a a X \ ^mod(exp) (Qa = ) 

i mod(exp)lV ; a — yv ) 

where 7 mo( j and 7 exp refer to the detector output spectrum calculated through the model de- 
scribed by Eqs. Q-Q, and obtained directly from the experiment, respectively. From Eqs. 
([SJ-Q, it is clear that the p -parameter is a function of the sample ellipsometric parameters 
OF, A), azimuthal angle of the polarizer <j> p , and the extinction ratio and phase retardance of 
both polarizer and analyzer. It is worth noting that under ideal conditions, i.e. t\p(a) — and 
Pj = [1/ y/2 1 / y/2], and for (f) p =45°, the ellipsometric parameters of the sample can be di- 
rectly extracted from the p -parameter as 

^ = tan- 1 |p|, A = arg(p), (6) 
where arg{.) denotes the phase operator. 

4. Measurement and Calibration Procedure 

In order to accurately measure the ellipsometric parameters OF, A) and consequently the optical 
constants through THz-TDSE, one needs to know the polarization state vector of the detector 
Pj, extinction ratio and phase retardance of the polarizer/analyzer as well as their azimuthal 
offset angles A(j) p /A(j) a with respect to the /?-axis. Azimuthal offset angles may arise due to the 
uncertainty in the polarizer (analyzer), or the misalignment of the components (e.g. surface tilt 
of the sample), and can be easily modeled by replacing R((f) p ^) with R((f) p ^ +A(f) p ^) in Eq. 

The complex vector Pj can be obtained experimentally over the desired frequency range 
through a calibration scheme reported elsewhere l30l . Using a gold mirror as the sample in Fig. 
|4| and by adapting the calibration scheme in l30l for the reflection mode, it is easy to show that 



where 



S d = arg(K), (9) 

(10) 



/exp(^=0°)-/exp(^ = 90 o 



/exp(^=0°)+/exp(^ = 90 o ) 

Once Pj is known from Eq. ([7]) and through the reflection measurement from a gold mir- 
ror, the ellipsometric parameters are obtained through a regression algorithm which we have 
adapted from the conventional regression calibration method developed by B. Johs for rotat- 
ing element optical ellipsometers [31], and used here for the first time for THz-TDSE. In our 
proposed regression algorithm, the p -parameter defined in Eq. [5] is obtained experimentally 
(Pexp) from an isotropic sample at one or several different azimuthal angles of the polarizer 
(<j> P i,i = 1,...,A/). Then, the p -parameter from the model (p m od) described by Eqs. ([2])-([5j> is 
fitted to that from the experiment over the desired frequency range and simultaneously for all 
the polarizer azimuthal angles, using the Levenberg-Marquardt nonlinear regression algorithm. 
We define an error function Err based on least squares fitting as 

^ = LL0Pmod(^^^,^,A,...)|-|p exp ((O i -,^)l] 2 + (11) 

(j) pi (Oj 

[^(Pmod(^^^,^,A,...))-arg(p exp ((O i -,^))] 2 - 

It should be noted that all measurements and computations in this work are made in discrete 
time and frequency domains. In the regression algorithm, the ellipsometric parameters \F and A, 
as well as the polarizer/analyzer parameters Tjp^, S P ^ and A(j> p ^ are simultaneously found 
such that the error function defined in Eq. ( pT| ) is minimized over the measurement frequency 
range (C0jj = 1, ...,M) and for all measured data sets corresponding to different polarizer an- 
gles ((j) p i, i = 1, ...,AT). The advantage of this technique is that the polarizer/analyzer parameters 
do not need to be a priori knowledge, although reasonable initial values are necessary for the 
convergence of the regression algorithm. It was found that the best initial values for *P and A 
are those obtained from Eq. ([6]) when p = p exp . 

5. Experimental results and discussion 

As a demonstration, two different silicon samples were tested. The first sample was a 
0.4-mm-thick high resistivity Si substrate with polished surface. For this sample the reflected 
time-domain pulse was truncated right before the arrival of the second reflection in order to 
avoid Fabry -perot interferences due to the reflection from the back surface. The second sample 
was a 1-mm-thick highly phosphorus-doped Si substrate with nominal room temeprature DC 
resistivity of 0.0151 £2-cm. 

Figure [5] shows the measured calibrated and uncalibrated ellipsometric parameters *F and A 
for the high-resitivity sample. The calibrated plots are obtained by fitting the p -parameter to 
three sets of measured data corresponding to polarizer angles (j) p \ = 25°, (f) p 2 = 45° and (j) p 3 = 
65°. The uncalibrated plot is directly obtained from Eq. ^ using the data set corresponding 
to (f) p 2 = 45°. An incidence angle of 45° (away from the Brewster's angle) was used for this 
experiment so as to have appreciable signal in both the s— and /?— polarizations. 

Figure [6] compares the extracted refractive index, using Eq. ([I]), from the uncalibrated and 
calibrated ellipsometric parameters shown in Fig. [5] with that from the conventional transmis- 
sion THz-TDS reported in Q. It is evident that the calibration has substantially improved the 
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Fig. 5. Measured ellipsometric parameters (a) *F and (b) A for the high resistivity silicon 
substrate. Square and circle marks present the ellipsometric parameters before and after 
calibrating out the effect of the THz detector and the polarizer/analyzer, respectively. The 
incidence angle was 45°. 



measurement accuracy over the uncalibrated data. The ellipsometric measurements do not have 
as much precision as the transmission experiments. But this illustrates the important point 
that the power of ellipsometric measurements rests not on measuring dielectric properties of 
samples that could be measured in transmission, but in measuring samples that are otherwise 
opaque. 



-THz-TDS 

' THz-TDS E w/ calibration 
THz-TDSE w/o calibration 



0.7 0.9 
Freq. (THz) 



Fig. 6. Extracted refractive index of the high resistivity silicon substrate from the uncali- 
brated (square) and calibrated (circle) ellipsometric parameters. Solid line shows the data 
from the conventional transmission THz-TDS reported in (T]. 



Figure [7] shows the simultaneous fitting results of the amplitude and phase of the p mo d to p exp 
for the highly doped silicon sample over the frequency range 0.1-1.5 THz. The fitting was per- 
formed simultaneously for three sets of data corresponding to three different azimuthal angle of 
the polarizer. Ellipsometric parameters \F and A along with polarizer/analyzer extinction ratios, 
phase retardances and offset angles were the free parameters of the regression algorithm in this 
least squares fitting. This highly doped silicon sample is opaque in the THz range and cannot 



be measured in transmission. The incidence angle was 73°, which is close to the Brewster angle 
of silicon. 
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Fig. 7. Fitting results of the (a) amplitude and (b) phase of the p -parameter shown as solid 
lines for p mo( j and discrete marks for p exp . The fitting was performed simultaneously for 
three sets of data corresponding to three different azimuthal angle of the polarizer and over 
the displayed frequency range. The sample was a highly doped silicon substrate, and the 
incidence angle was 73°. 

Figure[8ja) shows the extracted complex refractive index of the highly doped silicon substrate 
from ellipsometric parameters (not shown) along with a Drude model fit to the measured data. 
Fig. [SJb) shows the AC resistivity of the highly doped silicon substrate extracted from the 
ellipsometry measurement. As shown, the measured AC resistivity is consistent with the DC 
resistivity measured via a non-contact eddy-current resistivity gauge (COTS ADE 6035). 

In order to avoid water vapor absorption, the space with terahertz wave propagation was 
enclosed and purged with dry air during the above measurements. 

6. Conclusion 

Terahertz time-domain spectroscopic ellipsometry is a newly established characterization tech- 
nique that collectively possesses the advantages of the well established techniques of ellip- 
sometry and THz-TDS. In this work, a new experimental setup for THz-TDSE was presented 
that easily provides variable angle of incidence in the range of 15°-85° in the reflection mode. 
The setup can also be configured into transmission mode with the same ease. We used hollow 
core photonic band gap fiber without pre-chirping to deliver femtosecond laser to the THz pho- 
toconductive antenna detector that facilitates the change of the incidence angle. An effective 
calibration scheme was proposed and applied for characterization of two different Si substrate 
samples as an initial demonstration. 
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Fig. 8. (a) Extracted complex refractive index of the highly doped silicon substrate from 
ellipsometric parameters. Solid lines show the Drude model fit to the measured data, (b) 
Extracted resistivity of the highly doped silicon substrate from ellipsometry measurement. 
Square mark presents its DC resistivity measured via a non-contact eddy-current resistivity 
gauge (COTS ADE 6035). 



